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Interacting fluids generating identical, dual and phantom cosmologies 
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We find the group of symmetry transformations generated by interacting fluids in spatially flat 
Friedmann-Robertson-Walker (FRW) spacetime which links cosmologies with the same scale factor 
(identity) or with scale factors a and a -1 (duality). There exists a duality between contracting and 
superaccelerated expanding scenarios associated with (phantom) cosmologies. We investigate the 
action of this symmetry group on self-interacting minimally(conformally) coupled quintessence and 
fc-essence cosmologies. 


I. INTRODUCTION 

Symmetry transformations preserving the form of the 
spatially flat FRW equations introduce an alternative 
concept of equivalence between different physical prob¬ 
lems pj. This means that a set of cosmological models 
are equivalent when their dynamical equations are form 
invariant under the action of that group. It suggests that 
any of this equivalent cosmologies can be used to de¬ 
scribe the present accelerated expansion of the universe. 
Hence, it turns out to be interesting to investigate the 
consequences of this group. 

Due to the additivity of the stress-energy tensor a cos¬ 
mological model with one fluid in flat FRW spacetime 
can be seen as a model of two interacting fluid compo¬ 
nents. For instance a scalar field minimally coupled to 
gravity has been described as a stiff fluid interacting with 
vacuum energy 0 ] . The stress-tensor of the tachyon field 
could be considered as the sum of two components, one 
behaving like a pressureless dust and the other having a 
negative pressure j 2 ]. In this letter we investigate and ex¬ 
tend these results splitting the source in two components 
in a manner compatible with two discrete symmetries of 
the Einstein equations in flat FRW spacetimes. These 
symmetries are the structural invariance of the scale fac¬ 
tor as a function of the cosmological time, i.e., a[t) = a(t) 
0,0,l,i and the duality between expanding and con¬ 
tracting backgrounds a = 1/ a 0,1,0 M- I n addition, 
the dual transformation mapping contracting into super- 
accelerated expanding backgrounds provides the link be¬ 
tween a standard and a phantom cosmology 1,0- A 
phantom source with sufficiently negative pressure vio¬ 
lates the weak energy condition but could describe ade¬ 
quately current observations [h|. Other characteristics 
of phantom cosmologies have been investigated in Ell¬ 
in section II we develop the interacting framework 
and illustrate it with simple examples, such as, self¬ 
interacting minimally(conformally) coupled quintessence 
field and fc-essence field <p. In section III we present 
a linear transformations which preserve the form of the 
dynamical equations and apply them to those scalar 


fields. In section IV the conclusion are stated. 


II. INTERACTING FRAMEWORK 

The Einstein equations in the flat FRW spacetime with 
scale factor a and a perfect fluid read 

3 H 2 = p, ( 1 ) 

p + 3 H(p+p)= 0 , ( 2 ) 

where p is the energy density, p the pressure and H = 
a/a. There are two independent Einstein equations for 
the three quantities a,p, and p. Usually, the system of 
equations ifTO) is closed with an equation of state p = 

p(p)- 

We assume that Tik splits into two perfect fluid parts, 
Tik = T ik + T iki with T lk = (Pl,2 +Pl, 2 )UiU k + Pipgik , 
where pip, and pq2 are the energy density and the equi¬ 
librium pressure of fluids 1 and 2 respectively, u l is the 
four-velocity. Therefore, Eqs. CEP become 

3 H 2 = pi + p- 2 , ( 3 ) 

Pi+ P2 + 3 H (pi + p2 + Pi + P2) = 0 . ( 4 ) 

The Eq. P shows the interaction between the fluid com¬ 
ponents allowing the mutual exchange of energy and mo¬ 
mentum. Consequently, there will be no local energy- 
momentum conservation for the fluids separately. To 
preserve the degree of freedom of the original system 
of equations up, we introduce an equation of state 
for each fluid component pq2 = (71,2 — l)pq2, where 
71,2 are the barotropic indexes of fluids 1 and 2 respec¬ 
tively. As the energy-momentum tensor of the system as 
a whole is conserved, we assume an effective perfect fluid 
description with equation of state p = (7 — l)p, where 
7 = {I1P1+I2P2) / P is the effective barotropic index. For 
this effective perfect fluid the dynamical equations are 
identical to CP - 

The scalar field </> with energy density and pressure 

P<l> = ^ 2 + v ( 4 >)^ P<t> = ~ V { 4 >), ( 5 ) 


can be represented in terms of two interacting fluids 
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Pi = p\ and P2 = ~P2, meaning that 71 = 2 (stiff mat¬ 
ter) and 72 = 0 (vacuum energy). Due to the interactions 
between the two fluid components the energy-momentum 
tensor conservation of the system as a whole is equivalent 
to the Klein-Gordon equation 


For the fc-essence field ip with Lagrangian C = 
— U(ip) F(x), where U(ip) is the potential, F(x) is a func¬ 
tion of the kinetic term x = g lk piPk and ipi = dip/dx *, 
we associate the energy-momentum tensor of a perfect 
fluid. The energy density and pressure of the k field are 


<j) + 3 Ft(p + — 0, (6) 

while 7 = 2pi/p0 is the effective barotropic index. 

The conformal scalar field ip driven by the potential 

V(V0 = + Vo, Vo > 0, (7) 


is an interesting model, because this potential has re¬ 
ceived much attention in the literature in connection with 
the early inflationary epoch 0- This simplified model 
leads to a final accelerated expansion phase retaining the 
essentials of minimally coupled approaches. The energy 
density and the pressure of the conformal scalar field 

Pp = ^(ip + Hip) 2 + Xip 4 + V 0 , ( 8 ) 

Pip = ^(.'•P + Hip ) 2 + ^ip 4 - Vo, ( 9 ) 

can be represented as the sum of two interacting fluids 
namely pi = (ip + Hip) 2 / 2 + A ip 4 and P2 = Vo, with 
equations of state pi = pi /3 and P2 = —P2 representing 
(radiation) 7! = 4/3 and (vacuum energy) 72 = 0 , while 
7 = 4 pi/ 3 pp is the effective barotropic index. The energy 
conservation of the effective fluid is given by the Klein- 
Gordon equation 

ip + 3Hip + -Rip + 4Xip 3 = 0 , ( 10 ) 
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whose first integral is 

^(ip + Hip) 2 + Xip 4 = \, ( 11 ) 

2 a 4 

where b is an integration constant. After combining Eqs. 
0 , 0 , JTTJ and integrating, we get the scale factors 


a 


± 


C 


a c 



b> 0 


b< 0. 


( 12 ) 

(13) 


Defining the conformal time 77 = f dt/a and the new field 
X = ipa, the Eq. m becomes x" + 4 Ay 3 = 0, where 
' = d/dip Its general solution can be expressed in terms 
of the Jacobi functions. For A > 0 , the qualitative aspects 
of ip are obtained from the first integral CD. which reads 
X ' /2 /2 + Ay 4 = b. Hence, x oscillates between — 6 1 / 4 < 
X < b 1 ! 4 and ip = x/a becomes a decreasing oscillating 
function with a vanishing final limit at the minimum of 
the potential. On the other hand, the bouncing solution 
m avoids the initial singularity. 


p v = U(tp)[F - 2 xF x ], Pv = —U(ip)F, ( 14 ) 

with F x = dF/dx. They can be split as two interact¬ 
ing perfect fluids such that pi = U(p)F(x) and p 2 = 
— 2 U(ip)xF x with equations of state pi = —pi andp2 = 0 . 
They play the role of (vacuum energy) 71 = 0 and (dust) 
72 = 1 . The energy-momentum tensor conservation of 
the effective fluid is the fc-field equation 

U' 

[F x + 2 xF xx ]f> + 3 HF x tp + — [F - 2 xF x ] = 0 , ( 15 ) 

while 7 = P2Ip<p is the effective barotropic index. 


III. FORM INVARIANCE SYMMETRY 

Here, we will find a symmetry transformation that pre¬ 
serves the form of the system of equations 0 and 0 . 
To begin with, we observe that the total energy density 
is form invariant under the linear transformations 

(£)-(i“ 0 V)(£). <*» 

that is, p = pi p2 = pi + P2 = P for any a and f 3 . 
These form invariant transformations constitute a group 
and they induce the transformations H = H or H = —Ft 
in Eq. 0 . The former leads to the identity a = a and 
the latter to the duality a = 1 /a. The duality between 
contracting and superaccelerated expanding ( H > 0 and 
Ft > 0 ) phases gives rise to a phantom transformation 
[5J. Assuming equations of state for the interacting flu¬ 
ids, pip = (71,2 — l)pi,2 in the initial configuration and 
Pi,2 — (71,2 — l)(h,2 in the final configuration, then the 
conservation equation 0 remains form invariant when 

72-71 a 72-71 ^ zj 

a,. = -—, Pi=~ -—, H = H , ( 17 ) 

72 -71 72 - 7i 

ad= f 2 ±Ji^ = H = -H. ( 18 ) 

72 -7i 72 - 7i 

The duality Ft —> —H, acting on Eq. 0 induces the 
transformation p + p —> — (p+p) and the associated mat¬ 
ter violates the weak energy condition. Such model could 
explain current observations and it is referred to as phan¬ 
tom cosmology [HI. When both p and p diverge the dual 
transformation trades a final big crunch by a final big 
rip. However, when p is finite but p diverges 0 , the 
dual transformation interchanges this finite time singu¬ 
larity by other of the same kind. These sudden future 
singularities can occur even when the matter obeys p > 0 
and p + 2 >p > 0. 
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Below, we show the action of the form invariance sym¬ 
metry group on the three fields <j>, 1 / and p. The as¬ 
sociated effective perfect fluid description will be used 
splitting each fluid as two interacting ones with indexes 
(71,72), i- e., = (2,0), C(i/) = (4/3,0) and 

k(p) = (0,1). Also, we allow real or imaginary fields. 

k —> Q: Let us suppose that we are interested in seek¬ 
ing the set of k-essence models having the same scale 
factor that the quintessence ones. To this end we use 
Eqs. ©, itTfll . (t!T)ll and 03- Then, we get 

= p 2 UF\ V = U[F - F'x\ . (19) 

We are demanding that the last be an identity, i.e. it does 
not impose any relation between the fields </, p and their 
derivatives. Under this requirement the second equation 
(HU translates into conditions 

V(<t>) = U(<p), F — F'x = 1. (20) 


The phantom sector of the duality comprises the (—) 
branch of Eq. CT for Q —> k and the branch ( 1251 ) for 
k —> Q. The former represents a cosmology with a fu¬ 
ture big rip singularity at t = 0 and a real(imaginary) p 
according to m > 0(to < 0). The latter describes a non 
singular cosmology with H > 0 and an imaginary <j>. 

C —> Q: We are going to find the relationships between 
a conformal scalar field model driven by the potential © 
and a scalar field one under an identical transformation. 
Using Eqs. ©, iiHiuii . on and 03, we obtain 

\<t> 2 = ^(/P + Hi/) 2 -I- ^AV> 4 , (26) 

V{</) = \^ + Hi/) 2 + \\i/ 4 FV 0l (27) 

the potential V(cj>) = Vo + </> 2 /4 and the scale factors ill 21 
03- Besides, the potential V can be reconstructed as a 
function of </ and 


The first equation indicates that both potentials are the 
same when written as functions of cosmological time, 
but different when written as functions of the individual 
fields. Integrating the second Eq. ©3) we get F = 1+mx 
with m an arbitrary integration constant. Inserting F 
into the first Eq. 03 we find the following relationship 
between the fields <f> and p ID? 


4> = V2m j VUdp. 


( 21 ) 


It shows the identity between quintessence and k-essence 
models in flat FRW. For instance, the inverse square po¬ 
tential U oc p~ 2 correspond to an exponential potential 
V oc e-v^/v^Vo |D| and U oc p 2n to V oc 1 ). 

For k —> Q the dual transformation QH , 03 gives 


1 : r 


= -p 2 UF\ V = U[F -3 F'x\. 


( 22 ) 


Requiring that these equations be identities, we have 
V{<j>) = U{p) and F = 1 — mx 1 / 3 . Then, the first Eq. 
gives the transformation rule for the kinetic terms 


2mU 


-p 


2/3 


(23) 


For constant potentials V = U = Vo > 0, the first inte¬ 
grals of Eqs. © and 03 can be written as = 2b/a 6 
and p 2 / 3 = 2 ba 6 where the integration constants b and b 
transform as b = 2mVo&/3 (see Eq. Q3 )- After solving 
the Eq. ©, we express the duality as 


a = 


±\j sinh y/SVot 

\j cosh y/SVot 


-1/3 


1 - 1/3 


1 

a = -, 
a 


1 

a = -, 
a 


b > 0, 

(24) 

b < 0. 

(25) 



2 2. 9 

cosh d> -sinh“ 6 
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correspond to the solutions 01 and 


(28) 


= -2iarctaneV' 4V °/ 3t , V = V 0 


cos 2 i(f>+ - sin 2 iq 

O 


(29) 


to the solutions 03- 
For C —> Q the dual transformation 03, 03 gives 


+ Hi /) 2 - ^A?/ 4 , 


(30) 


V(0) = ^(iP + Hi/) 2 + ^Xi/ A + V 0 . (31) 

These equations lead to the potential V((j>) = Vo — 5</> 2 /4, 
so comparing the first integral of the Eq. © with Eqs. 
03 and QH , we have c/ 2 = 2ba‘ i with b = —26/3. After 
solving the Eq. 0, we express the duality as 


at = {a±) 1 , b < 0, 

(32) 

a c = (a c ) _1 , b > 0, 

(33) 


The phantom sector of the duality comprises the (—) 
branch of Eq. 03 for C —► Q and the branch QH for 
Q —> C. The former represents a cosmology with a future 
big rip singularity at t = 0, where </> and V are obtained 
from II28I29I) making the substitution </> —> i(f>. The latter 
describes a non singular cosmology with U > 0. 

C —> k: From Eqs. (18191) . QH , and 0303 , we have 


UF = -hi/ + Hi/) 2 ~\\i/ A FV 0 , 
o 6 

-2 UxF’ = h/jj + Hi/) 2 + ^Xi/ 4 . 

o o 


(34) 

(35) 
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These equations lead to U = Vo and F = 1 + mi 2 , which 
along with Eqs. 0 and m give the general solution 
cast - For m < 0, the singular solution m of this 
purely kinematics k-essence model interpolates between 
radiation and de Sitter phases. For m > 0 the non sin¬ 
gular solution m bounces at t = 0. 

For G —» k the dual transformations m- ns gives 

UF=l^ + H^) 2 + h^ + V 0 , (36) 

-2 UxF' = + H VO 2 - ^\ip\ (37) 

with U = Vo and F = 1 + mx 2 ^ 7 . In this case the solu¬ 
tion of Eq. m becomes Eqs. (tiro and the duality is 
expressed by Eqs. m . 

C —> C: We complete this section investigating the du¬ 
ality between conformally coupled quintessences. Duality 
among Q —> Q or k —> k were studied in Refs. §j- Jj}|. 
From, Eqs. 00 , m and ©, we have 

^(ip + Hip) 2 + Xip 4 = -^{ip + Hipf - Xip 4 , (38) 

V 0 = V 0 + {ip + Hip) 2 + 2A ip A . (39) 

Comparing Eqs. (tTTI and (l3%l) we get 6/a 4 = —6/a 4 = 
6 a 4 , so 6 = 6 = 0 and V 0 = Vo by Eq. (BH1) . The Eq. 0 
turns into 3 H 2 = V 0 and a = e H ° t , with H 0 = ±^/V 0 /3. 
Integrating the Eq. we obtain the conformal fields 


ip = 


Hn 


1p = ~- 


Hn 


V= 2 A[l-e ff ot]’ r ^ 2 l[l 
along with their transformation rule 


iP = 




iP- 


Hn 


V^2A L 


(40) 

(41) 


So, the duality induces a linear transformation group act¬ 
ing on ip and there is no phantom sector because H = 0. 


IV. CONCLUSIONS 

We have found a symmetry group generated by the 
additivity of the stress-energy tensor and shown that the 
different forms of summing two interacting fluid compo¬ 
nents induces two discrete symmetries in the Einstein 
equations. They correspond to cosmologies with iden¬ 
tical geometry or to the duality between expanding and 
contracting backgrounds. The duality between contract¬ 
ing and superaccelerated expanding scenarios connects 
standard and phantom cosmologies. 


The identical transformation relates two different cos¬ 
mologies, both with the same scale factor, so the choice 
of a particular model to describe the evolution of the uni¬ 
verse is not unique because there are a variety of ways of 
combining two interacting fluids with the same geome¬ 
try. We highlight that in spatially flat FRW cosmologies, 
the dual transformation gives rise to linear relations be¬ 
tween the components of the Einstein tensor, which read 
Gq = Gq, G 4 = 2Gq — G\ and so on for the remain¬ 
ing spatial components. Consequently, the conservation 
laws preserve their form i.e., ViG^. = V,G/ = 0. This 
analysis reveals that no matter how accurately future ex¬ 
periments may come to determine the cosmological ob¬ 
servables there will always be a fundamental uncertainty 
about which of the possible models leading to the ob¬ 
served set of values has been chosen by Nature. This 
uncertainty could be avoided by considering other kind 
of frameworks as superstring field theory, M-theory or 
some specific interaction between the parts in which the 
total fluid was divided. For instance, in Ref. 0 it was 
shown that a particular interaction between both fluids 
reduces the linear transformations 0J) to the identity 
restricting the form invariance symmetry. 

We have shown that fc-essence cosmologies generated 
by a linear kinetic function and quintessence ones share 
the same scale factor. 

The conformally coupled quintessence model we have 
studied is equivalent to a purely kinematics k-essence 
model generated by a quadratic kinetic term. It describes 
a universe that begins to evolve as if it were radiation 
dominated at early times and ends in a constant vacuum 
energy dominated phase. 

For Q —> k, C —> Q and C —> k, we have found a 
phantom symmetry between contracting universes end¬ 
ing in a big crunch and expanding universes ending in a 
big rip. In these examples, the transformation rule for 
the scalar fields we have obtained is different than Wick 
rotation. In the case G —> G, the dual transformation is 
so restrictive allowing only the duality between expand¬ 
ing and contracting de Sitter geometries and inducing a 
linear relation among the conformal scalar fields. 

Finally, the linear transformations we have found can 
be extended to the case of fluids with non constant 
barotropic indexes and we think this is an interesting 
subject to be investigated in the future. 
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